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1.  INTRODUCTION 

Many  U.S.  Department  of  Defense  (DoD) 
installations  contain  soil,  sediment,  surface  water,  and 
groundwater  environments  contaminated  with 
explosives  and  energetics  (Brannon  and  Myers,  1997). 
Both  TNT  and  RDX  are  widely  used  explosives. 
Human  exposure  can  occur  by  drinking  contaminated 
water,  breathing  contaminated  air,  or  coming  in  contact 
with  contaminated  soil.  Risk  assessments  should 
adequately  link  a  contaminant  source,  transmission 
pathways,  and  exposure  potential.  Meanwhile, 
compliance  with  water  quality  regulations  is  becoming 
increasingly  important  at  military  installations. 
Ecosystem  management  to  provide  for  sustained  and 
future  sustainable  mission  capacity  remains  the  DoD 
policy  for  military  installations.  The  ecosystem 
management  expectation  can  only  be  met  by  providing 
the  tools  necessary  to  actively  manage  watersheds. 
Water  quality  and  related  aquatic  ecosystems  are  major 
endpoints  and  are  insufficiently  understood  components 
of  natural  resource  management,  particularly  for 
military  installations.  These  installations  should  be 
provided  with  immediately  usable  and  effective  models 
that  can  be  implemented  for  compliance  with  water 
quality  regulations  as  well  as  long-term  watershed 
planning  and  management.  Modeling  as  part  of 
watershed  management  to  meet  water  quality  goals  is 
not  new  but  most  current  models  were  developed  and 
tested  two  decades  ago.  Watershed  models  are  largely 
confined  to  lumped  and  semi-distributed  surface  water 
simulation.  Models  that  reflect  hydrologic  and  aquatic 
impacts  from  military  conditions  are  rare.  Watershed 
models  that  enable  diagnostic,  predictive,  and 
operational  applications  in  conjunction  with  monitoring 
and  data  collection  programs  are  virtually  non-existent 
across  the  board  and  are  urgently  needed  within  the 
scientific  and  modeling  communities. 

To  meet  this  need  a  physically  based,  distributed 
source  Contaminant  Transport,  Transformation  and 
Fate  (CTT&F)  sub-module  was  developed  by  the  U.S. 
Army  Engineer  Research  and  Development  Center,  to 


simulate  point  and  non-point  sources  across  a 
watershed.  The  model  operates  on  a  cell  by  cell  basis. 
Soil  types,  land  uses,  and  other  hydrologic 
characteristics  can  be  varied  spatially  among  cells. 
Potential  chemicals  are  routed  through  the  cells  from 
the  watershed  divide  to  the  outlet.  The  distributed, 
process-oriented  structure  of  the  model  can  be  used  for 
identifying  critical  source  areas  of  non-point  sources 
within  the  watershed,  and  can  give  insight  on  both  the 
persistence  and  fate  of  explosives.  The  expectation  is 
that  a  model  of  this  type  is  able  to  quantify  transport 
and  transformation  of  multiple  contaminants  and  can 
facilitate  the  assessment  of  the  fate  of  distributed 
sources  and  lead  to  better  management  of  the  watershed 
environment  for  military  installations.  CTT&F  can  be 
linked  to  any  distributed  hydrologic  model,  assuming 
the  hydrologic  model  provides  the  required  flow  and 
sediment  transport  fluxes.  The  ability  of  the  model  to 
simulate  explosives  in  watersheds  is  demonstrated  by  a 
test  case  study  in  the  laboratory  environment. 

2.  GOVERNING  EQUATIONS 

To  simulate  the  chemical  transport  processes  in 
watersheds,  it  is  necessary  to  estimate  beforehand  the 
hydrological  variables.  The  hydrological  variables 
required  to  drive  CTT&F  can  be  calculated  using  any 
distributed  model  capable  of  producing  a  reasonable 
simulation  of  the  flow  and  sediment  transport  fields  in 
the  watershed.  As  presented  in  the  SHETRAN  model 
(Ewen,  et  ah,  2000),  the  basic  model  components 
include  hydrology,  erosion,  and  sediment  and  chemical 
transport.  The  calculations  for  each  process  at  any  time 
level  are  independent  and  information  is  carried 
forward  from  the  hydrology  sub-module  to  the 
sediment  transport  sub-module  to  the  chemical 
transport  sub -module  in  order  to  generate  a  solution. 

2.1  Hydrology  and  flow  sub-module 

Considerable  advances  have  been  made  in 
physically  based  distributed  hydrologic  modeling 


1 


Report  Documentation  Page 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 

1 .  REPORT  DATE  2.  REPORT  TYPE 

01  NOV  2006  N/A 

3.  DATES  COVERED 

4.  TITLE  AND  SUBTITLE 

Development  of  a  Distributed  Watershed  Contaminant  Transport, 
Transformation  and  Fate  (CTT&F)  Sub-Module  for  Military 

Installations 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Environmental  Laboratory,  U.S.  Army  Engineer  Research  and 
Development  Center,  CEERD-EP-W  Vicksburg,  MS  39180 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release,  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

See  also  ADM002075.,  The  original  document  contains  color  images. 

14.  ABSTRACT 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF:  17.  LIMITATION  OF 

18.  NUMBER  19a.  NAME  OF 

a.  REPORT  b.  ABSTRACT  c.  THIS  PAGE 

unclassified  unclassified  unclassified 

38 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


watershed  in  the  past  few  years.  There  are  several  grid- 
based,  distributed  models  to  simulate  watershed 
hydrologic  and  sediment  transport  processes.  The  U.S. 
Army  Corps  of  Engineer’s  GSSHA  (Gridded  Surface 
Subsurface  Hydrologic  Analysis)  is  a  physically-based, 
distributed-parameter,  structured  grid,  hydrologic 
model  that  simulates  the  hydrologic  response  and 
sediment  transport  of  a  watershed  subject  to  given 
hydrometeorological  inputs  (Downer  and  Ogden, 
2004).  GSSHA  is  a  reformulation  and  enhancement  of 
the  CASC2D  model  and  encompasses  the  full 
hydrologic  cycle. 


governed  by  conservation  of  mass  (sediment  continuity) 
(Julien,  1998): 


dc,  |  dqu  |  (k!n 
dt  dx  dy 


+  W, 


(3) 


where  Cs  =  concentration  of  sediment  particles  in  the 
flow  [M/L3];  qa,  qlx  =  total  sediment  transport  areal  flux 
in  the  x-  or  y-direction  [M/L2T];  Je  =  sediment  erosion 
volumetric  flux  [M/L3T];  Jd  =  sediment  deposition 
volumetric  flux  [M/L3T];  and  Ws  =  sediment  point 
source/sink  volumetric  flux  [M/L3T]. 


The  governing  equations  for  water  flow  with  the 
GSSHA  or  similar  distributed  models  are  based 
primarily  on  the  Saint- Venant  equations  of  continuity 
and  momentum.  Overland  flow  is  routed  in  two 
orthogonal  directions  (x  and  y)  in  each  grid  cell  during 
each  time  step.  Once  water  enters  a  "channel"  grid  cell, 
then  the  volume  of  water  is  added  to  the  channel  system 
and  routed  to  the  watershed  outlet.  The  2D  continuity 
equation  for  gradually-varied  flow  over  a  plane  is: 


dh  dq  v  dq 

—  +— —  +  — - 
dt  dx  dy 


(1) 


where  h  =  surface  water  depth  [L];  qx,  qy  =  unit 
discharges  in  the  x-  or  y-direction  [L2/T];  and  ie  = 
excess  precipitation  rate  [L/T] . 


The  diffusive  wave  approximation  of  the  Saint- 
Venant  equations  have  found  wide  application  in  many 
areas  and  are  well  established  for  modeling  a  variety  of 
watershed  hydrologic  processes.  The  2D  diffusive  wave 
momentum  equations  are  described  as: 


5  -  S  -  — 

fx~  0x  dx  ’ 


Sfy  ~  $0y 


dll 

dy 


(2) 


where  S/X,  Sfy  =  friction  slope  in  the  x-  or  y-direction; 
and  Sfh,  S(lx  =  ground  surface  slope  in  the  x-  or  y- 
direction. 


The  governing  equations  for  the  channel  flow 
routing  process  are  similar  to  those  for  overland  flow, 
the  flow  and  suspended  sediment  transport  in  channels 
is  described  by  the  1-D  equation. 

2.3  Chemical  transport,  transformation  and  fate 
sub-module 

The  CTT&F  sub-module  includes  all  possible 
physical  and  biochemical  processes  taking  place  in  the 
watershed.  Physical  transport  mechanisms  affect  the 
spatial  variability  of  chemical  mass  throughout  the 
watershed.  Partitioning  and  biochemical  reactions, 
meanwhile,  determine  the  distribution  of  chemical  mass 
among  different  phases  and  thereby  affect  the  amount 
of  mass  available  for  transport.  In  consideration  of 
these  important  processes,  the  governing  equations  for 
chemicals  are  established  over  a  differential  control 
volume  through  which  the  fluid  is  flowing.  A  basic 
principle  of  contaminant  transport  models  is  the 
conservation  of  mass.  Mathematically  modeling  for 
chemical  transport  involves  the  simultaneous  solution 
of  governing  equations  of  both  water  flow  and  upper 
sediment.  The  soil  sediment  bed  plays  an  important  role 
in  the  transport  of  contaminants.  Particles  and 
associated  chemicals  in  the  surficial  sediments  may 
enter  deeper  sediment  layers  by  burial  or  be  returned  to 
the  water  columns  by  scour.  An  overview  of  the 
processes  included  in  CTT&F  is  presented  in  Fig.  1. 


2.2  Erosion  and  sediment  transport  sub-module 

Sediment  erosion  and  transport  processes  take 
place  both  on  the  land  surface  and  within  the  channel 
network.  Sediment  erosion  and  transport  are  potentially 
very  important  processes  at  the  watershed  scale.  Excess 
sediment  can  affect  water  quality  directly.  Sediment 
transport  also  influences  chemical  transport  and  fate. 
Many  chemicals  sorb  strongly  to  sediment  and  thus 
undergo  settling,  scour,  and  sedimentation.  Sorption 
also  affects  a  chemical's  mass  transfer  and 
transformation  rates.  Both  sediment  transport  rates  and 
concentrations  must  be  simulated  in  most  toxic 
chemical  studies.  For  the  overland  plane  in  2D,  the 
concentration  of  sediment  particles  in  a  flow  is 


Fig.  1.  Schematic  chart  of  the  key  processes  of  CTT&F 
sub-module 


2 


2.3.1  Chemical  partitioning  and  distribution 

CTT&F  incorporates  “four  phase  partitioning”, 
where  a  contaminant  may  partition  between  solid 
particles,  sediments,  DOC,  and  water.  Pure  solid 
chemicals  are  assumed  to  be  a  separate  phase  and  are 
modeled  as  reactive  solid  particles.  Equilibrium 
partitioning  between  the  truly  dissolved  (aqueous), 
particle-sorbed,  and  DOC-bound  phases  is  assumed. 
The  equilibrium  partition  coefficient  can  be  used  to 
describe  the  fraction  of  the  total  chemical  that  is 
associated  with  each  phase.  The  concentration  of  the 
chemical  in  truly  dissolved,  bound,  and  particulate 
phases  then  can  be  calculated  by: 


Cc  =  fdCT ,  Cb  =  fbCT,  Cp  =  ZfpncT  (4) 

f  = _ kbCDoc _  (5a) 

1  +kbCDOC  +  ZkpnCpn 

f  = _ kpnCpn _  (5b) 

P"  1  +  kbCDOC+YJkpnCpn 

/,=i-A-Z4„  (5c> 

where  Cd  =  dissolved  phase  chemical  concentration 
[M/L3];  Cb  =  bound  phase  chemical  concentration 
[M/L3];  Cp  =  particulate  phase  chemical  concentration 
[M/L3];  CT  =  total  chemical  concentration  of  the 
dissolved,  DOC  bound  and  sorbed  phases  [M/L3];/rf  = 
fraction  of  the  total  chemical  in  the  dissolved  phase;  fb 
=  fraction  of  the  total  chemical  in  the  bound  phase;  fp  = 
fraction  of  the  total  chemical  in  the  particulate  phase 
associated  with  particle  “n”;  kb  =  equilibrium  binding 
coefficient  [L3/M];  CDOC  =  concentration  of  DOC 
[M/L3];  kpn  =  particle -dependent  partition  coefficient 
[L3/M];  and  Cpn  =  concentration  of  particle  “n”  [M/L3]. 

2.3.2  Chemical  Transport 

Overland  transport  of  chemicals  is  vital  for 
quantifying  a  distributed  source.  The  transport 
processes  that  are  included  in  the  overland  flow  in 
CTT&F  are  advection,  dispersion,  infiltration,  erosion 
and  deposition,  and  mass  transfer  between  overland 
flow  and  the  upper  soil.  Consideration  of  these 
important  processes  leads  to  the  2D  advection- 
dispersion  equation  (ADE).  The  governing  equations 
for  the  total  chemical  concentration  in  the  overland 
flow  and  upper  soil  layer,  respectively,  can  be  written 
as  follows  (Johnson  and  Zhang,  2006). 


Overland  flow: 

d(hCrT)  |  d{qxCrT)  |  d( qyCrT )  d 
dt  dx  dy  dx 


hD 


dc; 

dx 


(6) 
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hD 


dc; 
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-  ~J  df  +J  dm+J  pe  ~  J  nd  +  Z  5 


pd 


Upper  soil  layer: 

d(~CT2)  _  .  _  .  _  .  .  ,y  c  (7) 

g  J  df  J  dm  Jpe+Jpd+^-‘‘->k  V  ' 

where  CrT  =  total  chemical  concentration  in  the 
overland  flow  [M/L3];  C'T2  =  total  chemical 
concentration  in  the  upper  soil  layer  [M/L3];  z  =  depth 
of  the  upper  soil  layer  [L];  D„  D,  =  chemical  dispersion 
coefficient  in  the  x-  or  y-direction  [L2/T];  Jd;  = 
dissolved  chemical  infiltration  flux  [M/L2/T];  Jdm  = 
dissolved  chemical  mass  transfer  flux  between  surface 
water  and  upper  soil  layer  [M/L2/T];  Jpe  =  sorbed 
chemical  erosion  flux  [M/L2/T];  Jpd  =  sorbed  chemical 
deposition  flux  [M/L2/T];  and  E Sk  =  total  biochemical 
transformation  flux,  positive  indicates  a  source, 
negative  a  sink  [M/L3/T], 

Similarly,  two  conceptual  areas  are  defined  within 
the  channel:  the  main  channel  and  the  bed  sediment. 
Internal  transport  processes  of  a  chemical  occur  in  the 
channel  are  advective  and  dispersive  movement, 
sedimentation,  resuspension,  and  mass  transfer  between 
bulk  water  and  benthic  sediments.  Consideration  of 
these  important  mechanisms  leads  to  the  familiar  ID 
ADE.  The  governing  equations  for  the  total  chemical 
concentration  in  the  channel  flow  and  benthic  sediment, 
respectively,  can  be  written  as  follows  (Johnson  and 
Zhang,  2006). 


Channel  flow: 

dc;  +  QX  dc; 


dt 

<h 


A  dx 


]_d_ 
A  dx 


AD, 


dc; 

dx 


(8) 


(c;-c;)+Ujp 


pd 


■  J  dt )  +  J  dd  +^Sk 


Bed  sediments: 

d(Zf/2)=-Jpe  +  JPd-Jdd+ZSk  (9) 

where  C*  =  total  chemical  concentration  in  the 

channel  flow  [M/L3];  C.”2  =  total  chemical 

concentration  in  the  bed  sediments  [M/L3];  Qx  =  flow  in 
the  channel  direction  [L3/T];  J dt=  dissolved  chemical 

transmission  loss  flux  [M/L2T];  and  J dd=  pore  water 
diffusion  flux  of  dissolved  chemicals  [M/L3/T], 

Beyond  partitioning  and  mass  transport  processes, 
transformation  of  explosives  has  been  observed  as 
summarized  by  Townsend  and  Myers  (1996)  and  their 
fate  is  potentially  influenced  by  a  number  of 
biochemical  transformation  processes  that  include,  but 
are  not  limited  to,  biodegradation,  hydrolysis,  oxidation, 
photolysis,  volatilization.  These  processes  may  not  be 
operative  in  all  environmental  settings.  Some  chemicals 
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(12) 


undergo  a  complex  set  of  reactions,  while  others 
behave  in  a  more  simplified  manner.  The  importance  of 
any  particular  process  is  highly  dependent  on  the 
chemical  of  interest  and  environmental  settings.  The 
chemicals  may  be  independent  or  they  may  be  linked 
with  reaction  yields,  such  as  a  parent 
compound-daughter  product  sequence.  CTT&F  allows 
the  simulation  of  a  variety  of  biochemical  processes 
that  may  affect  chemicals,  which  were  adapted,  in  part, 
from  the  WASP  (Ambrose,  et  al.,  1993)  and  IPX 
(Velleux,  et  al.,  2000)  water  quality  models  because 
biochemical  transformation  mechanisms  and  their 
interrelations  are  poorly  understood  for  explosive 
compounds.  The  transport  equations  described  above 
must  be  supplemented  by  the  following  biochemical 
process  fluxes  as  source  or  sink  terms. 

2.3.3  Biodegradation  of  chemicals 

CTT&F  simulates  biodegradation  as  a  first-order 
kinetics  process,  which  means  that  the  rate  of  loss  of 
mass  at  any  given  time  is  directly  proportional  to  the 
mass  present  at  that  time,  the  rate  coefficients  can  be 
combined  into  a  single  decay  coefficient. 

3%-  =  kt,.Cr  <10 

at 

where  kbio  =  overall  biodegradation  rate  [1/T]. 

2.3.4  Hydrolysis  of  Chemicals 

Hydrolysis  is  a  reaction  in  which  cleavage  of  a 
molecular  bond  of  the  chemical  and  formation  of  a  new 
bond  with  either  the  hydrogen  or  the  hydroxyl 
component  of  a  water  molecule  occurs.  Hydrolysis  by 
specific  acid,  neutral,  or  base-catalyzed  pathways  is 
considered  for  the  various  species  and  phases  of  each 
chemical.  The  total  rate  of  hydrolysis  transformation  of 
a  chemical  is  computed  by  CTT&F  as  the  sum  of  three 
contributing  processes. 

khyd=Tkjln+]+Ki+kbJ[0H- ])/.  (ID 

j 

where  khyd  =  overall  hydrolysis  rate  [1/T];  kaj,  kbj  = 
specific  acid  catalyzed  and  base  rate  constants  for 
chemical  in  phase  j,  respectively  [L3/M/T];  knj  =  neutral 
rate  constant  for  chemical  in  phase  j  [1/T];  and  fj  = 
fraction  of  chemical  in  phase  j. 

2.3.5  Oxidation  of  Chemicals 

Oxidation  is  modeled  by  CTT&F  as  a  general 
second-order  process  for  the  various  species  and  phases 
of  each  chemical.  Oxidation  reactions  are  computed 
from  the  chemical  input  data  and  the  total 
environmental  concentrations  of  reactive  oxidizing 
species. 


kox,=\R02\LkOJfJ 

i 

where  koxi  =  overall  oxidation  (or  reduction)  rate 
[1/T];  [RO2]  =  molar  concentration  of  oxidant  [M/L3]; 
and  k0j  =  second  order  oxidation  rate  constant  for 
chemical  in  phase  j  [L3/M/T]. 

2.3.6  Photolysis  of  chemicals 

Photolysis  or  photodegradation  is  the 
transformation  of  a  chemical  due  to  absorption  of  light 
energy.  The  first  order  rate  coefficient  for  photolysis  is 
calculated  by  CTT&F  from  the  absorption  rate  and  the 
quantum  yield  for  chemical  in  each  phase. 

kphl=l.kjjfj  (13) 

j 

where  kpi„=  overall  photolysis  rate  [1/T];  kaj  =  specific 
sunlight  absorption  rate  for  chemical  in  phase  j  [E/M/T]; 

and  (/)  ■  =  reaction  quantum  yield  for  chemical  in  phase  j 

[M/E], 

2.3.7  User-defined  extra  reaction  of  chemicals 

An  extra  user-defined  second-order  reaction  for  the 
various  phases  of  each  chemical  is  included  in  CTT&F. 
The  second  order  reaction  allows  the  user  to  simulate 
the  effect  of  processes  not  considered  by  the  model. 
The  reaction  depends  upon  a  rate  constant  and  an 
environmental  parameter. 

Krd=[E]Zkejfj  (14) 

i 

where  kurrj  =  net  user-defined  reaction  rate  constant 
[1/T];  [£]  =  intensity  of  environmental  property  driving 
the  reaction;  and  kei  =  second  order  rate  constant  for 
chemical  in  phase  j,  [E]"  T'  . 

2.3.8  Dissolution  of  solid  phase  chemicals 

Explosives  associated  compounds  are  commonly 
present  as  crystalline  solids  (McGrath  (1995). 
Dissolution  of  solid  phase  chemicals  into  water  is  a 
primary  mechanism  by  which  solid  chemical 
contamination  spreads  through.  Once  dissolved,  the 
chemicals  are  available  for  transport  or  other 
biochemical  transformation  processes.  For  explosive 
chemical  compounds  dissolution  rates  have  not  been 
widely  studied.  In  CTT&F,  the  first-order 
approximation  of  the  dissolution  process  is  used  to 
describe  explosive  dissolution  rate  (Cussler,  1997): 

=kdAS-(fd+fb)CT]  (15) 
at 

where  Cs  =  solid  phase  chemical  concentration 
[M/L3];  kdsl  =  the  dissolution  mass  transfer  coefficient 
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[L/T] ;  a  =  area  available  for  mass  transfer  between  the 
solid  and  liquid  [L2];  and  S  =  aqueous  solubility  of  the 
chemical  [M/L3]. 

2.3.9  Transformations  and  daughter  products 

The  chemicals  simulated  by  CTT&F  may  be 
independent,  or  they  may  be  linked  with  reaction 
yields,  such  as  a  parent  compound-daughter  product 
sequence.  Linked  transformations  may  be  implemented 
by  simulating  two  or  three  chemicals  and  by  specifying 
appropriate  yield  coefficients  for  each  process. 

(16) 

j  k 

where  Jkji  =  production  of  chemical  "i"  from  chemical 
"j"  undergoing  reaction  "k"  [M/L3/T];  Kkj  =  effective 
rate  coefficient  for  chemical  "j"  reaction  "k"  [1/T];  and 
Ykji  =  yield  coefficients  for  production  of  chemical  "i" 
from  chemical  "j"  undergoing  reaction  "k"  [M/M]. 

2.3.10  Numerical  solutions 

The  coupled  set  of  CTT&F  differential  equations  is 
solved  by  numerical  techniques,  which  uses  a  finite 
difference  (FD)  control  volume  solution  scheme.  A 
watershed  system  is  discretized  into  a  mesh  of  square 
grids,  the  locations  of  which  are  described  in  terms  of 
rows,  columns,  and  layers.  DEM-derived  local  drainage 
directions  are  used  as  the  basis  for  channel  routing. 
However,  the  channel  routing  is  conducted  up  and 
down  slope  lines  across  the  watershed.  By  using  a  FD 
algorithm,  solution  of  the  model  yields  a  general 
equation  of  the  form: 

C(x,  v, t )  =  ( Hydrologic  Transport)+  (Transformation  )  (17) 


3.  MODEL  TESTING  AND  VALIDATION 
3.1  Experiment  set-up 

Given  required  data  from  field  studies,  CTT&F  can 
be  used  to  assess  the  fate  of  explosive  compounds  in 
watershed  systems.  Efforts  to  model  RDX  and  TNT 
from  Camp  Shelby  fire  range  watershed,  Mississippi, 
the  largest  state-owned  training  site  in  the  nation,  are 
currently  underway  with  the  CTT&F  sub-module. 
Unfortunately,  reliable  watershed  field  data  is,  up  to 
now,  unavailable.  In  order  to  validate  the  general 
performance  of  the  model,  the  CTT&F  sub-module  was 
evaluated  by  means  of  a  laboratory  test  case  study  of 
explosives  transport  and  transformation  in  a  laboratory. 
A  laboratory  experiment  was  conducted  to  simulate 
rainfall  and  overland  flow  associated  with  sediment  and 
chemical  transport.  The  experimental  procedure  was 
designed  to  mimic  field  conditions  resulting  from  the 


direct  action  of  rainwater  on  target  distributed  sources 
in  military  installations.  The  experiment  plot  was  9.0’x 
7.5’.  The  bed  slope  of  the  plot  was  2%  and  designed  to 
collect  runoff  water  and  sediment  from  the  plot.  Two 
land  covers  were  used  as  part  of  the  experiments  to 
simulate  two  different  surface  roughness:  “disturbed” 
(unvegetated)  and  “undisturbed”  (vegetated).  The  soils 
used  in  the  experiment  were  obtained  from  the  Camp 
Shelby  military  fire  range.  The  chemical  and  physical 
properties  of  the  soils  before  rainfall  were  analyzed  by 
standard  test  methods.  The  rainfall  was  introduced 
through  a  rainfall  simulator.  Rainfall  intensity  for  the 
overall  plot  area  was  ranged  from  2.7  to  2.9  in/hr.  The 
simulated  rainfall  event  lasted  30  +  60  +  90  min. 
Runoff  and  sediment  were  collected  at  the  downstream 
end  of  the  plot  and  measured  volumetrically.  Runoff 
rates  and  volumes  were  determined  by  collecting  runoff 
every  minute  of  a  30  minute  rainfall  simulation  and 
every  minute  after  rainfall  was  discontinued  until  it  was 
noted  that  runoff  had  ceased.  The  total  suspended 
sediment  (TSS)  samples  were  collected  every  minute 
for  the  first  15  minutes  of  runoff,  then  every  five 
minutes  during  the  30  minutes  rainfall  simulations  and 
every  minute  afterward.  To  simulate  the  transport  and 
transformation  of  RDX  and  TNT,  a  total  of  500  g 
Comp  B  of  varying  sizes  was  applied  uniformly  onto 
the  soil  surface.  After  the  Comp  B  application,  the  soil 
surface  was  subjected  to  a  simulated  rainfall  event, 
which  induced  overland  flow  and  chemical  transport. 
Comp  B,  which  has  been  widely  used  in  munitions,  is 
commonly  present  as  crystalline  solids  and  is  a  60/39 
mixture  of  RDX  and  TNT  that  contains  1%  wax.  TNT 
solubility  is  greater  than  RDX.  The  chemical  reaction 
and  transport  caused  by  each  rainfall  event  were 
measured  by  collecting  samples.  During  each  rainfall 
event,  4-  liter  runoff  samples  were  collected  every  5 
minutes  for  chemical  analysis  of  explosives;  their 
concentrations  were  determined.  This  test  case 
primarily  explores  the  dissolution  and  overland 
transport  capabilities  of  CTT&F. 

The  experiment  plot  was  modeled  using  a  domain 
consisting  of  30  grid  cells,  the  size  of  the  grid  was  1 .5  ft 
by  1.5  ft.  The  time  interval  for  the  simulation  was  0.5 
sec.  The  physical,  chemical,  and  biological  properties 
of  RDX  and  TNT  and  their  kinetic  reaction  rates 
required  by  the  model  were  estimated  by  using  the  data 
collected  from  the  experiment.  Table  1  shows  the 
actual  average  explosive  compound  concentrations  for 
three  Comp  B  samples.  Other  values  for  the  relevant 
parameters,  as  showed  in  Table  2  were  adopted  from 
Lever  et  al.  (2005)  and  McGrath  (1995),  but  ultimately 
they  were  adjusted  by  trial-and-error  to  reproduce  the 
measured  concentrations  of  RDX  and  TNT  from  the 
experiment.  Many  transformation  processes  affecting 
explosives  are  active  primarily  in  the  aqueous  phase 
and  therefore  may  be  limited  by  dissolution  kinetics. 
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This  study  gives  the  results  in  predicted  solid 
persistence  and  multi-phase  chemical  loading  from 
solids  into  aqueous  solution  for  TNT  and  RDX. 

Table  1.  Analysis  for  three  Comp  B  particles 
Comp  B  HMX  RDX  TNT 

(mg/kg)  (mg/kg) (mg/kg) 

1  59424  562798  350955 

2  68039  637121  393580 

3  71505  672170  422214 


Table  2.  Model 

parameters 

for  RDX  and  TNT  in  water 

Parameter 

Comp 

B 

RDX 

TNT 

Aqueous  solubility 

- 

4.6  x  10'5,a 

1.3  xl0'4,a 

(g/cm3) 

Diffusion 

_ 

2.2  x  10'6'a 

6.7  x  10'6'a 

coefficient(cm2/s) 
Density  (g/cm3) 

1.65a 

1.82b 

1.654b 

1st  transformation 

" 

0-  1.0  x  10"1,b 

- 

rate  (1/hr) 

“from  Lever  et  al.  (2005) 
bfrom  McGrath  (1995). 


3.2  Results  and  discussion 

Numerical  results  were  obtained  from  running  the 
CTT&F  sub-module.  In  this  experiment,  overland  flow 
causes  erosion  and  dissolution  of  the  solid  Comp  B,  a 
fraction  of  which  infiltrates  into  the  soil  while  the 
remainder  is  transported  downstream.  Even  though 
distributed  observations  for  RDX  and  TNT 
concentrations  were  not  measured  in  this  study,  we  can 
infer  and  trace  the  migration  of  distributed  RDX  and 
TNT  sources  using  the  model. 

In  this  study,  the  model  was  calibrated  by 
comparing  the  simulated  and  measured  surface  runoff, 
sediment  concentration,  and  chemical  concentration 
such  that  the  errors  between  simulated  and  measured 
values  were  minimum.  The  calibration  results  for  this 
test  case  study  using  the  selected  model  parameters  are 
listed  in  Tables  4.  Methods  used  to  evaluate  model 
performance  included  comparison  of  simulated  and 
measured  results  through  statistical  analysis.  The  most 
fundamental  approach  to  assessing  model  performance 
in  terms  of  behavior  is  through  visual  inspection  of  the 
simulated  and  observed  results.  To  describe  the  model 
goodness  of  fit  to  the  overall  hydrographs  and 
pollutographs,  the  RMSEs  (Root  Mean  Square  Errors) 
were  calculated  and  given  in  Tables  4.  The  R2  and 
Nash-Sutcliffe  Efficiency  (NSE)  criteria  for  assessing 
the  efficiency  of  the  model  are  also  given  in  Tables  4. 
For  surface  runoff  from  both  unvegetated  and  vegetated 
plots,  the  simulated  values  compared  reasonably  well 
with  measurements.  The  event  averaged  percent  errors 


of  both  simulated  total  surface  discharges  were  less 
than  10%  of  its  corresponding  measured  value.  The 
RMSE  and  R2  values  between  simulated  and  measured 
results  for  the  unvegetated  plot  were  1.195  and  0.723, 
respectively.  For  the  vegetated  plot  the  RMSE  and  R2 
values  between  simulated  and  measured  results  were 
0.641  and  0.944,  respectively.  The  model  did  not  fully 
capturing  the  initial  wash-off  of  sediments,  the  event 
averaged  percent  error  of  simulated  TSS  concentration 
from  both  unvegetated  and  vegetated  plots  was  31.76% 
and  55.22%,  respectively.  The  RMSE  was  considered 
to  be  high  and  the  R2  value  was  low.  The  model 
performance  for  suspended  sediment  concentration  was 
strongly  affected  by  the  initial  six  samples  collected, 
and  the  extremely  high  sediment  concentrations  were 
measured  from  these  samples.  Furthermore,  the  model 
was  capable  of  capturing  the  general  trends  of  TSS 
concentration  over  time,  the  simulated  values  compared 
well  with  measured  values  after  the  initiation  of  the 
event.  The  Error,  RMSE,  R2,  and  NSE  values  are 
greatly  improved  without  the  inclusion  of  the  first  six 
samples.  For  explosive  compounds,  simulated  errors 
are  very  small  (within  7%).  The  R2  values  between 
simulated  and  measured  concentration  results  from  the 
unvegetated  plot  were  0.995  and  0.997  for  RDX  and 
TNT,  respectively.  The  R2  values  between  simulated 
and  measured  results  from  the  vegetated  plot  were 
0.687  and  0.895  for  RDX  and  TNT,  respectively. 
Further,  the  model  performed  well  consistently  for  two 
different  data  sets.  Comparisons  of  the  overall  shape  of 
simulated  and  measured  results  over  time  for  surface 
runoff  discharge,  TSS  concentration,  dissolved  RDX 
and  TNT  concentrations  in  surface  runoff  from  (a) 
un vegetated  plot  and  (b)  vegetated  plot  are  shown  in 
Fig.  2  to  5,  respectively.  As  can  be  observed,  the 
agreement  of  model  simulations  and  experimental 
results  for  the  explosive  compounds  from  the  field  is 
satisfactory  thus  showing  that  the  CTT&F  sub-module 
is  able  to  capture  the  essence  of  explosive  transport  and 
transformation  processes.  Although  not  investigated  in 
this  study,  the  behavior,  transport  and  ultimate  fate  of 
distributed  sources  in  watersheds  may  be  affected  by 
other  chemical  reactions  and  related  phenomena. 
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Fig.  2.  Comparison  of  simulated  and  experiment  results 
over  time  for  surface  runoff  discharge 


Fig.  4.  Comparison  of  simulated  and  experiment  results 
over  time  for  RDX  concentration  in  overland  flow 


Fig.  3.  Comparison  of  simulated  and  experiment  results 
over  time  for  TSS  concentration  in  overland  flow 


Fig.  5.  Comparison  of  simulated  and  experiment  results 
over  time  for  TNT  concentration  in  overland  flow 


4.  SUMMARY  AND  CONCLUSIONS 

A  contaminant  transport,  transformation  and  fate 
sub-module,  CTT&F,  has  been  developed  for 
simulating  distributed  sources  in  watersheds.  CTT&F 
is  able  to  simulate  both  surface  runoff  and  channel 
processes  of  contaminants.  The  CTT&F  equations  are 
comprehensive,  physically  based,  and  fully  compatible 
with  various  distributed  hydrologic  models  which 
provides  the  required  hydrological  and  sediment 
variables.  The  physical  basis  is  important  since  it 
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Table  4.  Comparison  of  simulated  and  measured  hydrographs  and  pollutographs 


Parameter 

Simulated 

Measured 

Error  (%) 

R2 

RMSE 

NSE 

unvegetated  plot 

Surface  runoff  (L/min) 

189.72 

201.75 

5.96 

0.723 

1.195 

0.685 

Total  suspended  sediment  (mg/L) 

20917.60 

30653.33 

31.76 

0.166 

719.47 

0.231 

Dissolved  RDX  (mg/L) 

2.805 

2.782 

0.84 

0.995 

0.012 

0.994 

Dissolved  TNT  (mg/L) 

3.806 

3.776 

0.79 

0.997 

0.012 

0.997 

vegetated  plot 

Surface  runoff  (L/min) 

151.20 

139.83 

8.13 

0.944 

0.641 

0.923 

Total  suspended  sediment  (mg/L) 

726.02 

2106.67 

65.53 

0.04 

134.00 

0.247 

Dissolved  RDX  (mg/L) 

1.155 

1.207 

4.32 

0.687 

0.052 

0.532 

Dissolved  TNT  (mg/L) 

0.443 

0.417 

6.34 

0.895 

0.014 

0.865 

provides  the  link  between  the  simulations  and  physical 
property  measurements..  The  model  computes  on  a  grid 
basis  for  considering  spatially  varied  soils,  land  uses, 
and  other  hydrologic  characteristics.  CTT&F  can 
quantify  the  changes  in  environmental  conditions  such 
as  chemical  loads,  corresponding  to  various 
management  scenarios,  through  physically  meaningful 
parameters.  CTT&F  was  tested  to  demonstrate  its 
performance  using  an  experiment  plot.  Comparisons 
between  simulated  and  measured  results  have  been 
described.  The  comparisons  showed  that  the  model  was 
capable  of  simulating  the  explosive  compounds  from 
the  field  with  reasonable  accuracy.  Overall 
comparisons  were  encouraging,  and  showed  promise 
for  the  potential  use  of  the  CTT&F  sub -module  for 
studying  the  fate  of  distributed  sources  in  watersheds. 
Further  research  and  development  will  help  to  identify 
the  dominant  biochemical  transformation  processes, 
followed  by  better  descriptions  for  the  dominant 
processes  and  then  better  estimates  of  their  parameter 
values  for  each  variable  supporting  the  model.  The 
applicability  of  the  model  will  be  further  validated  and 
evaluated  using  more  data  sets. 
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Background 


Many  U.S.  Military  installations  potentially  contain  soil,  sediment,  surface  water, 
and  groundwater  environments  contaminated  with  explosives  and  energetics. 


The  release  of  contaminants  from  military  installations,  transport  across  the 
land  surface,  and  delivery  to  stream  networks  and  ground  water  can  have 
adverse  water  quality  and  ecological  impacts.  Compliance  with  water  quality 
regulations  is  becoming  increasingly  important  at  military  installations  thus 
necessitating  the  development  of  usable  and  effective  models  for  compliance 
as  well  as  long-term  watershed  planning  and  management. 
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Current  CTT&F  Sub-model’s  Features 


•  Spatial  Grid  Discretization 

•  Four  Phases  Partitioning: 

-  dissolved 

-  bound  to  DOC  (Dissolved  Organic  Carbon) 

-  sorbed  to  sediment  particles 

-  separate  solid  particles 

•  2D  Overland  Flow  Transport  with  Upper  Soil  Layer 

•  1 D  Channel  Flow  Transport  with  Bed  Sediments 


•  Seven  Biochemical  Transformation  Processes: 
-  biodegradation 


dimoMlor)  w  f  m  ild  p  ha  se 


-  user-defined  extra  reaction 
\  tran eforrn a tlon s  and  daughter  produeti 
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Chemical  Partitioning 

Equilibrium  partitioning  of  contaminants  among 
dissolved  phase,  sediment  sorbed  phase,  and  DOC 
bound  phase. 
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Biochemical  Transformation  Processes 


Biodeqredation|is  the  process  by  which  organic  substances  are 
broken  down  by  other  living  organisms. 

CTT&F  simulates  biodegradation  as  a  first-order  kinetics  process. 


Hydrolysis:  is  a  reaction  or  process  in  which  a  molecule  is  split  into 
two  parts  by  reacting  with  a  molecule  of  water,  which  has  the  chemical 
formula  H2O..The  reactions  are  first  order  for  the  neutral  chemical  and 

or  basic  forms  of  the  ch 


The  lorsiJ  !rsiie  of  bydrolyojo  IrsinoforTnsirjon  or  si  oherniosil  jo  computed 
by  GTTb.F  sio  the  oum  of  three  contributing  processes;. 
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Biochemical  Transformation  Processes 


Oxidation:  reactions  include  all  chemical  processes  in  which  atoms 
have  their  oxidation  number  (oxidation  state)  changed.  Oxidation 
describes  the  loss  of  an  electron  by  a  molecule,  atom  or  ion. 


Oxidation  is  modeled  by  CTT&F  as  a  general  second-order  process 
for  the  various  species  and  phases  of  each  chemical. 


Photolysis:  is  a  chemical  reaction  in  which  a  chemical  compound  is 
broken  down  by  photons  (e.g  linhn 


The  f\mi  order  mis  eoefriGieini  for  pbotolyete  te  calculated  by  CTTSa 
from  file  absorption  irate  end  the  quantum  yield  for  ohemioaJ  fa  each 
phase  7  j  ,  r 
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Biochemical  Transformation  Processes 


User-Defined  Extra  Reaction:  An  extra  user-defined  second- 
order  reaction  for  the  various  phases  of  each  chemical  is  included  in 
CTT&F. 


K.  =[E]Y.K,fj 


Biochemical  Transformation  Processes 


Dissolution  of  Solid  Phase:  Explosives  associated  compounds 


are  commonly  present  as  crystalline  solids.  Dissolution  of  solid  phase 
contaminants  into  water  is  a  primary  mechanism  by  which  solid 
contamination  spreads  through  the  system.  Once  dissolved,  the 
contaminants  are  available  for  other  transport  or  transformation 
processes.  In  CTT&F,  the  first-order  kinetic  process  is  used  to 
describe  the  explosive  dissolution  rate: 


dC, 

dt 


=  kd,,a[S-(fJ+ft)CT] 


niiim 
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Model  Testing  and  Validation  Studies 

•  The  test  plot  was  set-up  as  9.0’x  7.5’.  The  bed  slope  of  the  plot  was 
designed  with  2%  for  collecting  runoff  water  from  the  surface  of  the 
soil.  The  soils  were  obtained  from  Mississippi  Camp  Shelby  military 
fire  range.  The  rainfall  (2.8  in/hr)  was  introduced  through  a  rainfall 
simulator. 


•  Comp  B  is  a  60/39  mixture  of 
RDX  and  TNT,  Comp  B  is 
commonly  present  as 
crystalline  solids  and  were 
simulated  by  uniformly 
spreading  over  the  land 
surface 
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Model  Testing  and  Validation  Studies 


Model  parameters  for  RDX  and  TNT  in  water 

Parameter 

Comp 

B 

RDX 

TNT 

Aqueous  solubility  (25°C)  ( g/cm 3) 
Diffusion  coefficient  (25°C) 

(cm2/s) 

Density  (g/cm3) 

1st  order  transformation  rate  (1/hr) 

-1.65a- 

4.6  x  10 5a 

2.2  x  10 6  a 
1.82b 

0  -  1.0  x  10  lb 

1.3  xlO 4  a 
6.7  x  10  6  a 
1.654b 

a.  from  Lever  et  al.  (2005) 

b.  from  McGrath  (1995). 


Lever,  J.H.,  Taylor,  S.,  Perovich,  L.,  Bjella,  K.,  Packer,  B.,  2005.  “Dissolution  of  composition 
B  detonation  residuals.”  Environ.  Sci.  Technol.,  39(22):  8803-8811. 

McGrath,  C.  J.  1995.  “Review  of  formulations  for  processes  affecting  the  subsurface 
transport  of  explosives.”  Technical  Report  IRRP-95-2,  U.S.  Army  Engineer  Waterways 
Experiment  Station,  Vicksburg,  MS. 
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FY06  Model  Development  Activities 

Test  Plot  Model  Hydrology  Results 

Comparison  of  simulated  and  experiment  results  over  time  for  surface 
runoff  discharge. 

U  n  vegetated  Vegetated 


FY06  Model  Development  Activities 


Test  Plot  Model  Sediment  Results 


Comparison  of  simulated  and  experiment  results  over  time  for  total 
suspended  sediment  concentration  in  overland  flow. 

Un  vegetated  Vegetated 
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FY06  Model  Development  Activities 

Test  Plot  Model  Contaminant  Results 

Comparison  of  simulated  and  experiment  results  over  time  for 
dissolved  RDX  concentration  in  overland  flow. 


U  n  vegetated  Vegetated 
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FY06  Model  Development  Activities 

Test  Plot  Model  Contaminant  Results 

Comparison  of  simulated  and  experiment  results  over  time  for 
dissolved  TNT  concentration  in  overland  flow. 


U  n  vegetated  Vegetated 
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FY06  Model  Development  Activities 


Watershed  Validation  Site 


Camp  Shelby  Training  Sites 
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FY06  Model  Development  Activities 

Validation  Site  -  Camp  Shelby 

•  Camp  Shelby,  Mississippi,  the  largest  state  owned  training  site  in  the  nation, 
has  a  long  history  of  serving  the  country  and  is  considered  by  many  as  “a 
national  treasure.”  During  wartime,  the  camp's  mission  is  to  serve  as  a  major, 
independent  mobilization  station  of  the  U.S.  Army  Forces  Command 
(FORSCOM).  Camp  Shelby 


•  Training  Site  is  the  largest  reserve  component  training  site,  covering  136,000 


acres,  allowing  up  to  battalion  level  maneuver  training,  excellent  FA  Firing 
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Validation  Site  -  Camp  Shelby 
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FY06  Model  Development  Activities 


Camp  Shelby  Gage  Data 


Camp  Shelby 


Middle  Creek  (FY04) 
Stage  and  Flow  Data 
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FY06  Model  Development  Activities 


Camp  Shelby  Gage  Data 


Pearces  Creek  (FY05) 
Stage  Data 


SiDm  'I'evkuDiDglmfor  iha  'd'mmfomiml  Amsy** 


Camp  Shelby 
Firing 


Streams 

lmpact_area.shp 
Watershed  Boundary 
Cstsroadshp 

N 


USGS  02479144  PEARCES  CREEK  NR  JANICE,  MS 


0 

Mar  19  Mar  26  Apr  02  Apr  09 


nmental  Quality  and  Installations 


FY06  Model  Development  Activities 


Middle  Creek  Model  Development 
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FY06  Model  Development  Activities 


Middle  Creek-  Simulated  Hydrograph 
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FY06  Model  Development  Activities 


Middle  Creek-  Simulated  Distributed  Flow 
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FY06  Model  Development  Activities 


Middle  Creek-  Simulated  Contaminants  (RDX,  TNT) 


RDX 


TNT 


7.00E-12 


6.00E-12 


”  5.00E-12 


g  3.00E-12 

c 

o 

o 

X  2.00E-12 

Q 

CC 

1.00E-12 


1 

Total 

— ♦—  Dissolved 

1 

**•!+***+*+*  t« - X - T - 

7.00E-10 


6.00E-10 


«  5.00E-10 


o  3.00E-10 
o 
c 
o 

i-  2.00E-10 


1.00E-10 


Total 

— * —  Dissolved 

3/31/05  0:00  3/31/05  4/1/05  0:00  4/1/05  12:00  4/2/05  0:00  4/2/05  12:00 

3/31/05  0:00  3/31/05  4/1/05  0:00  4/1/05  12:00  4/2/05  0:00  4/2/05  12:00 
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FY06  Model  Development  Activities 

Middle  Creek-  Simulated  Distributed  Contaminants  -  RDX  and  TNT 
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FY06  Model  Development  Activities 


Model  Validation  -  Pearces  Creek 
Explosives  Flux  in  Surface  Water  During 
Storm  Events 

•  Water  quality  samples  collected 
during  two  storm  events  at  two 
gauged  streams  that  drain  central 
impact  area  at  Camp  Shelby 

•  Storms  represent  two  different 
discharge  conditions  -  enable 
calculation  of  explosives  flux 


^USGS 


USGS  02479144  PEARCES  CREEK  NR  JANICE,  MS 


^USGS 


USGS  02479144  PEARCES  CREEK  NR  JANICE,  MS 


nmental  Quality  and  Installations 


Summary  and  Conclusions 


•  A  CTT&F  sub-model  has  been  developed  based  on  the  mass 
conservative  form  of  the  coupled  system  of  2D  overland  flow  and  1 D 
channel  flow  for  watersheds.  The  CTT&F  equations  are  comprehensive, 
physically  based,  and  are  compatible  with  the  physically  based  (flow  and 
sediment),  distributed  watershed  hydrologic  models. 

•  CTT&F  works  on  a  grid  basis  for  considering  spatially  varied  soils,  land 
uses,  and  other  hydrologic  characteristics. 

•  CTT&F  generates  time  series  and  spatial  outputs  over  time. 


•  CTT&F  was  tested  for  RDX  and  TNT  from  a  test  plot  in  laboratory.  The 
comoarisons  showed  that  RDX  and  TNT  concentrations  in  overland  flow 


